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Abstract. The development of a dynamical model for investigating the nucleon 
resonances using the reactions of meson production from irN, ■yN, N{e, e'), and 
N{i/, I) reactions is reviewed. The results for the A (1232) state are summarized 
and discussed. The progress in investigating higher mass nucleon resonances is 
reported. 



1. Introduction 

The study of excited nucleon states (N*) has long been recognized as an important 
step towards developing a fundamental understanding of strong interactions. It is 
an important part of the effort to understand the structure of the nucleon since 
the dynamics governing the internal structure of composite particles, such as nuclei 
and baryons, is closely related to the structure of their excited states. Within 
the framework of Quantum Chromodynamics (QCD), a clear understanding of the 
spectrum and decay scheme of the N* states will reveal the role of confinement and 
chiral symmetry in the non-perturbative region. 

The N* states are unstable and couple strongly with the meson-baryon continuum 
states to form nucleon resonances in meson production reactions on the nucleon. 
Therefore the extraction of nucleon resonance parameters from the reaction data is 
one of the important tasks in hadron physics. By performing partial- wave analysis 
of pion-nucleon elastic scattering data mainly during the years around 1970, many 
7V*'s have been identified. From the resonance parameters listed by the Particle Data 
Group [1] (PDG), it is clear that only the low- lying N* states are well established while 
there are large uncertainties in identifying higher mass nucleon resonances. 

With the construction of high precision electron and photon beam facilities, the 
situation changed drastically in the 1990's. Experiments at Thomas Jefferson National 
Accelerator Facility (JLab), MIT-Bates, LEGS of Brookhaven National Laboratory, 
Mainz, Bonn, GRAAL of Grenoble, and Spring-8 of Japan have been providing new 
data on the electromagnetic production of tt, t], K, uj, (j), and 2tt final states. These 
data offer a new opportunity to to investigate N* properties, as reviewed in Refs.[2l[3]. 

In addition to analyzing the world's data of meson production from t:N, and 
A^(e, e') reactions, we need to interpret the extracted N* parameters in terms of QCD. 
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Figure 1. The total cross section data of meson production in 7p reaction. Left: 
1 - TT and 2- IT production are compared. Right: KY { K+A, K+TP , K°T.+ ), 
rip, and ujp production are compared with some of the 1 — tt and 2 — tt production 

There are two possibilities. The most fundamental way is to confront the extracted N* 
parameters directly with Lattice QCD calculations and QCD-based hadron structure 
models. Here the most challenging problem is to handle the contributions from the 
baryon continuum which are coupled with the reaction channels. The second one is 
to develop dynamical reaction models to analyze the meson production data. Here 
the reaction mechanisms and the internal structure of baryons are modelled by using 
guidances deduced from our understanding of QCD and many-year's study of hadron 
phenomenology. In this article, we give a review of the dynamical reaction models 
developed in Refs. gHHIllZllHlillinilllllllllllll]- Other approaches for investigating 
A^* states have been reviewed in Refs.[2j[3]. 

In practice, the dynamical reaction models describe the meson-baryon reaction 
mechanisms by using phenomenological Lagrangians which are constructed by using 
the symmetry properties, in particular the Chiral Symmetry, deduced from many- 
years' studies of meson-nucleon reactions. Starting from a set of phenomenological 
Lagrangians for mesons and baryons, one would ideally like to analyze the meson- 
baryon reaction data completely within the framework of relativistic quantum field 
theory. The Bethe-Salpeter (BS) equation has been taken historically as the starting 
point of such an ambitious approach. The complications involved in solving the BS 
equation in the simplest Ladder approximation have been known for long time. It 
contains serious singularities arising from the pinching of the integration over the 
time component. In addition to the two-body unitarity cut, it has a selected set of n- 
body unitarity cuts, as explained in great detail in Refs. [TH[Tn]. Thus it is extremely 
difficult, if not impossible, to apply the approach based on the Bethe-Salpeter equation 
to study N* states. 

Since 1990 the irN and reactions have been investigated mainly by using either 
the three-dimensional reductions [TJ of the Bethe-Salpeter equation or the unitary 
transformation methods [U [18]. These efforts were motivated mainly by the success 
of the meson-exchange models of NN scattering [19], and have yielded the meson- 
exchange models developed by Pearce and Jennings [201 j National Taiwan University- 
Argonne National Laboratory (NTU-ANL) collaboration [Hllll], Gross and Surva[^. 
Sato and Lee[4l[5], Juhch Group [24l[25l[26l [22], Fuda and his collaborators [l3|2l], and 
Utretch-Ohio collaboration [29l I30j . The focus of all of these dynamical models was on 
the analysis of the data in the A (1232) region. In this article, we will only review the 
model developed in Refs. [H [S] by using the unitary transformation method. We will 
also review its extension[Sl [7] to study the A (1232) excitation in neutrino-induced 
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N{v, Itt)N reactions. 

The main challenge of developing dynamical reaction models of meson production 
reactions in the higher mass N* region can be seen in Fig[TJ We see that two-pion 
photo-production cross sections shown in the left-hand-side become larger than the 
one-pion photo-production as the 7p invariant mass exceeds W ^ lA GeV. In the 
right-hand-side, KY ( fsT+A, K'^T,'^, rjp, and ojp production cross sections 

are a factor of about 10 weaker than the dominant tt^tt^p production. From the 
unitarity condition, we have for any single meson production process 7iV MB with 
MB = ttN, -qN, ioN, KA, KT, 

^{TMB.fN - T*^j^.ig) = ^ T^p B' ,MB PM' B'Tm' B' 
M'B' 

where pa denotes an appropriate phase space factor for the channel a. The large 
two-pion production cross sections seen in Figd] indicate that the second term in 
the right-hand-side of Eq.(IT]) is significant and hence the single meson production 
reactions above the A region must be influenced strongly by the coupling with the 
two-pion channels. Similarly, the two-pion production 7A'^ — > ttttN is also influenced 
by the transition to two-body MB channel 

i(r7r7rJV,7Af — r^Ar,7r7rAr) = ^ Tlj, g, .^^^ pM' B'Tm' B' ,'yN 

M'B' 

+ T*^pf .^^j^PT:TTNTTrTTN,jN ■ (2) 

Clearly, a sound dynamical reaction model must be able to describe the two pion 
production and to account for the above unitarity conditions. Such a model has 
been developed by using the unitary transformation method in Ref.0 and applied to 
investigate ttN elastic scattering [10], 7A'^ ttN reactions [Tl] ttN — > rjN reactions |1 2). 
and ttN ttttN reactions [13]. In this article, we will also review these results. 

This article is organized as follows. In section 2, we explain the unitary 
transformation method developed in Rcf. 3l"j using a simple model. The constructed 
model Hamiltonian for investigating N* states is given in section 3. The multi-channel 
multi-resonance reaction model developed in Refs.[4l [8] for calculating the meson- 
baryon reaction amplitudes is presented in section 4. In section 5, we give formula for 
defining the TV- TV* transition form factors and calculating the cross sections of pion 
production from ttTV, 77V, N{e,e'), and -/V(^/, /) reactions. The results in the A (1232) 
region and in the higher mass N* region are reviewed in section 6. A summary and 
discussions of future developments are given in section 7. 

2. Unitary Transformation Method 

The unitary transformation method was essentially based on the same idea of 
the Foldy-Wouthuysenth transformation developed in the study of electromagnetic 
interactions. It was first developed in 1950's by Fukuda, Sawada and Taketani [32], 
and independently by Okubo[33]. This approach, called the FST-Okubo method, 
has been very useful in investigating nuclear electromagnetic currents [SU [35] and 
relativistic descriptions of nuclear interactions [36l [3ll [38]. The advantage of this 
approach is that the resulting effective Hamiltonian is energy independent and can 
readily be used in nuclear many-body calculation. 
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To illustrate the unitary transformation method, we consider the simplest 
phenomenological Lagrangian density 

L{x) ^ Lo{x) + Li{x) , (3) 

where Lo{x) is the usual free Lagrangians with physical masses mjv for the nucleon 
field ipN and m,r for the pion field (jj^r, and 

Li{x) = i'N{x)T M,7TNi^N{x)4>T,{x) ■ (4) 

Here r^v TrAr denotes the physical ttNN coupling (~ f^rNN)- The Hamiltonian density 
H{x) can be derived from Eqs.(I3])-(|3])by using the standard method of canonical 
quantization. We then define the Hamiltonian as 



H = J H{x,t = 0)dx. (5) 

The resulting Hamiltonian can be written as 

H = Ho + Hi, (6) 

with 

Ho= j dk[EN{k)h\b^ + E^{k)a\a^l (7) 
Hi 



= / dkidhdkSik - h ~ fc2)[(rjv,.Ar(fci - k2)blb^^a^J + (h.c)], (8) 



where b^ and {b and a) are the creation (annihilation) operators for the nucleon and 
the pion, respectively. For simplicity, we drop the terms involving the anti-nucleon 
operator. Note that H along with the other constructed generators P, K, and J, 
as studied in Refs.[551 137] . define the instant-form relativistic quantum mechanical 
description of ttN scattering. We will work in the center of mass frame and hence 
the forms of these other generators of Lorentz group are not relevant in the following 
derivations. 

The essence of the unitary transformation method is to extract an effective 
Hamiltonian in a "few-body" space defined by an unitary operator U, such that 
the resulting scattering equations can be solved in practice. Instead of the original 
equation of motion H\a >= Ea\a >, we consider 

H'\a>=Ec,\a>, (9) 

where 

H' = UHU'' , (10) 

\a> =U\a> . (11) 

In the approach of Kobayashi, Sato and Ohtsubo'HT' (KSO), the first step is to 
decompose the interaction Hamiltonian Hj Eq.® into two parts 

Hi = Hf + Hf, (12) 

where Hf defines the process a ^ be with 771^ ^ 7ni) -f- ttIq which can take place in the 
free space, and Hf defines the virtual process with nia < mi, + me- For the simple 
interaction Hamiltonian Eq.®, it is clear that Hf — and Hf — Hj. 

The KSO method is to define an appropriate unitary transformation U to 
eliminate the virtual processes from transformed Hamiltonian H' . This can be done 
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systematically by using a perturbative expansion of U in powers of coupling constants. 
As a result the effects of 'virtual processes' are included in the effective operators in 
the transformed Hamiltonian. 

Defining U — exp{—iS) by a hermitian operator S and expanding U — l — iS+... , 
the transformed Hamiltonian can be written as 

H' = UHU^ 

= HQ + Hf + h'^ + [H(,,iS] + [Hi,iS] + ^^[[Hq,iS],iS] +•••. (13) 

To eliminate from Eq. (jl3p the virtual processes which are of first-order in the coupling 
constant, the KSO method imposes the condition that 

Hf + [Ha,iS]=Q. (14) 

Since i?o is a diagonal operator in Fock-space , Eq. (fT4|) clearly implies that iS must 
have the same operator structure of Hf and is first order in coupling constant. By 
using Eq. (fT4|) . Eq. lfTS]) can be written as 

H' = Ho + H'j, (15) 

with 

H'j = H[ + [Hf, iS] + i [Hf, iS] + higher order terms . (16) 

Since Hf , Hf , and S are all of the first order in the coupling constant, all processes 
included in the second and third terms of the H'j arc of the second order in coupling 
constants. 

We now turn to illustrating how the constructed Hj of Eq. pB]) can be used to 
describe the ttN scattering if the higher order terms are dropped. We consider the 
simple Hamiltonian defined by Eqs.®-® which gives Hf = and Hf = Ftv^ttAt- 
Our first task is to find S by solving Eg. lfTi)) within the Fock space spanned by the 
eigenstates of Hq 

Ho[N >= niNlN > , (17) 
Ho[k,p>= (E^ik) + ENip))\k,p> , (18) 
Ho[ki,k2,p>= {{E^{ki) + E^{k2) + EN{p))\kiM,P> , (19) 



For two eigenstates / and i of iJo, the solution of Eq. p4)) clearly is 

<;|,.S)|.>.-1|^. (20) 

Ef - Ei 

For the considered Hf — Tn^ttN we thus get the following non- vanishing matrix 
elements 

<k,mS)\N> = -F^,.^.(fc)— -4^|±^^ , (21) 

< N[{^S)[k',p' > = |^_±£)^r^ (fc') , (22) 

mN - E^(k') - En[p') 
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and 

i ^rc.^lP-^ -5{k'-k,)5{jf-k-p)TN,.Mik,) 

<M(.5)|a:„a:„,^> - ^ m f/T n'^'l'^V , , + (1^2) 

E-„{k) + £;Ar(p) - E.^(ki) - Ej,(k2) - -EAr(p) 

r (,„ •, -K^-ki)5{p -k2-p) , , , , 

With the above matrix elements and recaUing that Hf = and Hf = Vn^t^n for 
the considered simple case, the matrix element of the effective Hamiltonian Eq. (fT6|) 
in the center of mass frame (p = —k and = —k') is 

< k\H'j\k' > = i ^[(< k\TN^.N\I X I\{iS)\k' > 
I 

- <k\iiS)\I ><I\TN^^N\k' >]. (25) 

The only possible intermediate states are |/ >= |iV > +|7r(/ci)7r(fc2)A^(P/) >. By 
using Eas. (l^ - ([M)) we then obtain 

<k\H'j\k'> =v^'^{k,k')+v^''Hk,k')- (26) 

1 



where 



v^^\k,k') = -r%^^^{k)[ 



2 '^''^'"^ "E^{k)+ENik)-mN 
1 



E^{k') + ENik') - mN 



]r^,,^(fc') , (27) 



z;(")(^,fc') = ir^,.Ar(fc')[ ^ 



2 '^£;Ar(fc)-i;,(fc')-£;w(fc + fc') 

+ — ^ ]TN,.N{k) . (28) 

EN{k')~E^{k)~EN{k + k') 

Note that up to the same order Eg. ipS)) should have an additional term 
which is the one-pion-loop contribution to the single nucleon state. Such a mass 
renormalization term is dropped in practice, since it is part of the physical nucleon 
mass in the resulting effective Hamiltonian. If we treat this mass renormalization 
explicitly, we then will not get a solvable few-body problem, but a many-body problem 
which is as complicated as the original field theory problem. We also note that v^^^ of 
Eq. (|27|) is due to the intermediate "physical" nucleon state state |/ >= >. This is 
the consequence of the unitary transformation which eliminates the "virtual" nN N 
process. Here we see an important difference between u*^*' and the so-called nucleon- 
pole term from approaches based on some models based on the three-dimensional 
reduction of Bethe-Salpeter equations and the time-order perturbation theory [57]. 
There is no bare mass and energy-dependence in v^'^\ 

With the above derivations, the effective Hamiltonian Eq. (fT6| can be explicitly 
written as 

H' ^Ho + V, (29) 
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where 



iJo= j dk[EN{k)h\bj, + E^{k)a\af:], (30) 
V = f dkdk'[v^'\k,k')+v^''^k,k')]aibl^aj:,b^j:, . (31) 



To see the analytic properties of the reaction amplitudes based on the effective 
Hamiltonian Ea. (P5)) . let us first recall how the bound states and resonances are defined 
in a Hamiltonian formulation. In operator form the reaction amplitude is defined by 

^(^) = y+y^^^,^(^)^ (32) 

or 

(33) 

The analytic structure of scattering amplitude can be most transparently seen by using 
the spectral expansion of the Low equation Eq. p3|) 

V- < k'\V\<P,^ >< <^>e,\V\k > 



< k'\t{E)\k > = < k'\V\k > + 



^ E 



E-E' + ^e 

(34) 

where Ett is the threshold of the reaction channels, and \E'^'' are the discrete 
bound states and the scattering states, respectively. They form a complete set and 
satisfy 

> = > , (35) 

H'\¥+^ > = E'\¥+'> > . (36) 

Of course bound state energies are below the production threshold Eth- We now 
note that because of the two-body nature of V defined by Eq ([3T|) . Eq. ([35| has the 
one- nucleon solution H'\N >= Ho\N >= mi^\N >. But it does not contribute to the 
second term of Eq.([34l) because < t:N\V\N >= 0. Thus the amplitude Eq.([34l) does 
not have a nucleon pole which corresponds to bound state with a mass of physical 
nucleon and is formed by the physical N and tt of the starting Lagrangian Eq. ([3]). 
This is consistent with the experiment. Clearly, our approach is very different from 
the S-matrix approach which requires that the nN scattering amplitude must have a 
pole at i? = toat. Similar feature is also obtained by using the unitary transformation 
of Shebeko et al.[39l|40]. 

To end this section, we mention that the unitarity condition only requires that 
an acceptable model must have unitarity cut in physical region E > m^r + ttln- This 
is trivially satisfied in the the model defined by the effective Hamiltonian Eqs. ((30|l - 
(|3ip since the interaction V is energy independent. This is an important advantage 
in applying the method of unitary transformation to develop a multi-channels multi- 
resonances reaction models for investigating meson-nucleon reactions in the nucleon 
resonance region, as developed in Ref. [5] . In a model with an energy-dependent V such 
as the Julich model [57] the unitarity condition is much more difficult to satisfy, and 
the analytic continuation of the scattering t-matrix defined by Eqs. p4p to complex 
E'-plane is in general much more complex. 
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Figure 2. Basic mechanisms of the Model Hamiltonian defined in Eqs.(42)-(44). 



3. Model Hamiltonian 

With the unitary transformation method explained in section 2, it is straightforward 
to derive a model Hamiltonian for constructing a coupled-channel reaction model with 
7iV, -kN , TjN and mrN channels. Since significant parts of the ttttN production are 
known experimentally to be through the unstable states ttA, pN, and perhaps also 
aN, we will also include bare A, p and cr degrees of freedom in our formulation. 
Furthermore, we introduce bare N* states to represent the quark-core components of 
the nucleon resonances. The model is expected to be valid up to = 2 GeV below 
which three pion production is very weak. 

The starting point is a set of Lagrangians describing the interactions between 
mesons (M = 7, tt, rj, p,uj,a- ■■) and baryons {B — N, A, N* ■ ••). These Lagrangian are 
constrained by various well-established symmetry properties, such as the invariance 
under isospin, parity, and gauge transformation. The chiral symmetry is also 
implemented as much as we can. The considered Lagrangians are given in Ref.[H]. 
For completeness, we recall in Appendix A parts of these Lagrangians which were 
used in investigating the A (1232) resonance. 

By applying the standard canonical quantization, we obtain a Hamiltonian of the 
following form 



H = J h{x, t = 0)dx 

= Ho + Hi, (37) 

where h{x, t) is the Hamiltonian density constructed from the starting Lagrangians 
and the conjugate momentum field operators. In Eg.pTp. Hq is the free Hamiltonian 
and 

Hi = ^ ^MB^B' + ^ hM'M"^M , (38) 
M,B,B' M,M\M" 

where T mb^b' describes the absorption and emission of a meson(M) by a baryon(_B) 
such as TriV N and TriV ^ A, and hM' m" describes the vertex interactions 
between mesons such as tttt <-> p and 777 tt. 

Our main step is to derive from Eqs. ([57|) - ([551) an effective Hamiltonian which 
contains interactions involving TTTriV three-particle states. This is accomplished by 
applying the unitary transformation method up to the third order in interaction Hi 
of Eq. ([55)1 . The resulting effective Hamiltonian is of the following form 

H^ff^Ho + V, (39) 
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with 



Ho = Y.K^, (40) 



where Ka = ym-^ +P^^ is the free energy operator of particle a with a mass m^, 
and the interaction Hamiltonian is 

V = Tv+V22+v', (41) 



where 



Ty = {X^^X! ^N-^MB + Tn'-^tt-kn) + ^ ^Af-^TTTr} + {h.C.} , (42) 



N* MB W 



V22= ^ Vm'B',MB +V^7r ■ (43) 
MB.M'B' 

Here h.c. denotes the hermite conjugate of the terms on its left-hand-side. In the 
above equations, MB = 'yN,7TN,T]N,TTA, pN,aN represent the considered meson- 
baryon states. The resonance associated with the bare baryon state N* is induced by 
the vertex interactions Tn^^mb and Tn-^^t^ttN- Similarly, the bare meson states M* 
= p, a can develop into resonances through the vertex interaction hM'—tvTr- These 
vertex interactions are illustrated in Figllla). Note that the masses Af^, and to°^, of 
the bare states N* and A/* are the parameters of the model which will be determined 
by fitting the t:N and tttt scattering data. They differ from the empirically determined 
resonance positions by mass shifts which are due to the coupling of the bare states with 
the meson-baryon scattering states. It is thus reasonable to speculate that these bare 
masses can be identified with the mass spectrum predicted by the hadron structure 
calculations which do not account for the meson-baryon continuum scattering states, 
such as the calculations based on the constituent quark models which do not have 
meson-exchange quark-quark interactions. It is however much more difficult, but more 
interesting, to relate these bare masses to the current Lattice QCD calculations which 
can not account for the scattering states rigorously mainly because of the limitation 
of the lattice spacing. 

In Ea. (j43p . vm'B'.mb is the non-resonant meson-baryon interaction and Vt,.^ is 
the non-resonant tttt interaction. They are illustrated in Fig[2l^b). The third term in 
Eq. (|4ip describes the non- resonant interactions involving iniN states 

v' = V23 + W33 , (44) 

with 

V23 = ^[ivTr-KN,MB) + [h.c.)] , 
MB 

They are illustrated in Figl2l[c). All of these interactions are defined by the tree- 
diagrams generated from the considered Lagrangians. They are illustrated in Figl3]for 
two-body interactions vm' b'.mb and in Figl4]for v-^ttN.mb- In practice, we neglect Vt^tt 
and VTrjrN,-KTrN- We also only consider t;^^jv,7rAr and Vt^t^n,iN of VjrTTN,MB- These two 
interactions are illustrated in FiglH The calculations of the matrix elements of these 
interactions were explained in details in Ref . [8] . Here we only mention that the matrix 
elements of these interactions are calculated from the usual Feynman amplitudes with 
the energies of off-mass-shell particles in the intermediate states defined by the three 
momenta of the initial and final states, as specified by the unitary transformation 
methods. Thus they are independent of the collision energy E. 



Dynamical Models of the Excitations of Nucleon Resonances 



10 



B Bn B' 

vS v° 

Figure 3. Mechanisms for fM'S'.A/s of Eq. 1 143 I I: (a) direct s-channel, (b) 
crossed u-cliannel, (c) one-particle-excliange t-channel, (d) contact interactions. 



4. Multi-channels Multi-resonances Reaction Model 

Our next task is to derive a set of dynamical coupled-channel equations for describing 
jN,-!tN MB reactions within the model space N* ® MB © tt-kN. The starting 
point is the Lippman-Schwinger equation for the scattering T-matrix 

< a\T{E)\b >-< a\V\b > + < a\V l T{E)\b > , (45) 

E — Ho + le 



where the interaction V is defined from the effective Hamiltonian in Eqs. p9)) - ((44l) . 
We choose the normalization that the T-matrix is related to the S-matrix by 

< a\SiE)\b >= - ImSHPa - Pb) < a\T{E)\b > . (46) 

Since the interaction V, defined by Eqs. (|4T|) - (|44l) . is energy independent, it is rather 
straightforward to follow the formal scattering theory given in Ref. [41] to show that 
Eq. ([l5|) leads to the following unitarity condition 

< a\T{E) ~T^E)\b -27ri^ < a\T\E)\c > 6{E, - E) < c\T{E)\b>, (47) 

c 

where a, 6, c are the reaction channels in the considered energy region. 

We cast Eq. (1451) into a more convenient form for practical calculations. In the 
derivations, the unitarity condition Eq. (j47p must be maintained exactly. We achieve 
this rather complex task by applying the standard projection operator techniques [42j. 
similar to that employed in a study of ttNN scattering [l^. The details of our 
derivations are given in Appendix B of Ref. [5]. To explain our coupled-channel 
equations, it is sufficient to present the formula obtained from setting rjv*— ►ttttAt = in 
our derivations. Here we explain these equations and discuss their dynamical content. 

The resulting MB M'B' amplitude Tm'B'.mb in each partial wave consists 
of a non-resonant amplitude tM' B' .mb{E) and a resonant amplitude tfj,g, f^.jg{E) as 
illustrated in Figs. [n]and[Sl It can be written as 

Tm'B' ,Mb{E) = tM' B' ,Mb{E) + tMig, MB{E) . (48) 



NN NN AN AN N 

(a) (b) (c) (d) (e) 



NN NA AN AA N 

(f) (g) (h) (i) (i) (k) 

Figure 4. The considered fTrJV.TririV of D23 of Eq.(44). 
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t, 



. R 



t 



M'BMVIB 'M'B'.MB ^ M'B'.MB 

Figure 5. Graphical representations of Eq. 1 148 I I. 



t 



M'B',MB 



^M'BMVIB 



t 



M'B',MB 



Figure 6. Graphical representations of Eqs. I I49I I and I I57I I. 



Figure 7. Graphical representations of Eqs. I l52t - ll53l l. 



The second resonant term in the right-hand-side of Eq. is defined by 



'-M'B'MB 



(E)^ TN^^M'B'iE)[D{E)],,,TMB^N;iE), (49) 



N',N' 



with 



[D{E)-\,{E) = {E- M%.)5,,, - %,,{E) 



where M^, is the mass of a bare N* state, and the self-energies are 
^ij{E) — Y^^mb^n*{E)Gmb{E)Tnj^mb ■ 



MB 



(50) 



(51) 



In general, the bare states mix with each other through the off-diagonal matrix 
elements of the self-energies. The dressed vertex interactions in Eq. and Eq. 

((?T|) illustrated in Fig. [7] are (defining Tmb^n* = rjy.^^^^) 



^MB^N*{E) — TmB^N* + ^ ^M'B'^N*GM'B'{E)tM'B',MB{E) , 



AI'B' 



^N*^Mb{E) — rN*^MB + ^ tMB,M'B'{E)GM'B'{E)TN*^M'B' 



(52) 
(53) 



AI'B' 



The meson-baryon propagator Gai b in the above equations takes the following form 

1 



Gmb{E) 



E-Kb-Km- ^mb{E) + ie 



(54) 
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where the mass shift Yimb{E) depends on the considered MB channel. It is 
T,mb{E) = for the stable particle channels MB = ttN, r/N. For channels containing 
an unstable particle, such as MB = ttA, pN, aN, we have 

with 



T;Mb{E) = [< MB\gv— — ^, —g\,\MB >]un~ connected , (55) 

tj — iV^r — — Ajv + If- 



gV — Ta^U-JV + hp^TTH + ha^TTTT ■ (56) 

In Eq. ([55|) " un — connected" means that the stable particle, tt or N, of the MB state is 
a spectator in the ttttN propagation. Thus Timb{E) is just the mass renormalization 
of the unstable particle in the MB state. It is important to note that the resonant 
amplitude t^,j,g, mb{E) is influenced by the non-resonant amplitude tM'B'.MB{E), as 
seen in Eqs. ((49 | - ([53l) . 

The non-resonant amphtudes tM'B'.MB in Eg. ([48]) and Eqs. ([52 ]) - ([53 ]) are defined 
by the following coupled-channel equations 

tM' ,B' ,Mb{E) = Vm'B'.MB {E)+ J2 ' B' ,M" B" {E)Gm" B" {E)tM" B" ,M b{E) ,(57) 

M"B" 

with 

Vm'B' ,mb{E) — vm'B'.mb + Zm'B',mb{E) . (58) 

Here Zj^j>b\mb{E) contains the effects due to the coupling with ttttN states. It has 
the following form 



Zm'B',Mb{E) = [< M'B' I F- —Ft | MB >Unnected , (59) 

with 

(60) 

F = gv + vmb.tvttn , (61) 

where gv has been defined in Eq. (|55)l . Note that the dis-connected term in Eg. (155)1 
is already included in the mass shifts Sa/b of the propagator Eg. ([5^ and must be 
removed to avoid double counting. 

The appearance of the projection operator PttttN in Egs. (|55| and (|59|) is the 
consequence of the unitarity condition Eg. (|Tf)) . To isolate the effects entirely due to 
the vertex interaction gv = I'a^-kN + ^p^u-tt + ^o-^irTr , we use the operator relation 

^ + I (62) 

E- Ho-v E - Ho E-Hq E-Hq-v ^ ' 

to decompose the ttttN propagator of Eq. ((59)) to write 

Zm'B',Mb{E) = ^M's',Afs(-^) + ^AI'B'MBiE) ■ (63) 

The first term is 

Zi%',MBiE) = [< M'B' I gy—J^^f^gl \ MB >] connected • (64) 

Obviously, Zj^j/g, j^^.jgiE) is the one-particle-exchange interaction between unstable 
particle channels ttA, pN , and aN , as illustrated in FiglSj The second term of Eq. 
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z(E) = 

^'^'^^ N= N ^= A 



( E') ( E) ( E) 

Figure 8. One-particle-exchange interactions Z^^^^^{E), Z^^^^j^ and Z^^^^j^ 
ofEq. Jell. 



IS 



Z^^i^'B'MBiE) =< M'B' I F . t^^r,,.ME) ^ ^^T^^ . \ MB > 



hi — iio + ie 



+ < I VM'B'...N ^ ^z:"", . gi I MB > 

+ < M'B' I v,,,^,^,,^-I^^^^vl^ ,,,^ I MB > . (65) 
Here tT^7rN.7nTN{E) is a three-body scattering amplitude defined by 

(E) = W^^AT + V^-^N-^ -fp -r; Z —VttttN , (66) 

E - K„ - - A AT - Vj,T,N + «e 
where Vt^ttn has been defined in Eq. lpOl) . 

The amplitudes Tm'B',mb = ^m' b' ,MB+t^i' b' mb defined by Ea. (|15)) can be used 
directly to calculate the cross sections of -kN ttN, r/N and jN — > TriV, r/N reactions. 
They are also the input to the calculations of the two-pion production amplitudes. 
The two-pion production amplitudes resulted from our derivations are illustrated in 
FigM They can be cast exactly into the following form 

T-^7tN,mb (E) = T^I''n^mb (E) + T^:^pf jy.ig [E) + T^^jyjy.jg{E) + T^j^^./s {E) , (67) 
with 

'E^ZN,MBiE) =< V'Ljv(-^)I X! '"^^N.M'B'[Sm'B',MB 

M'B' 

+ GM'B'{E){tM'B',MBiE)+tfj,B,,MB)]\MB > , (68) 

T:XMBiE) =< V'i;]v(^)|ri^.^G.A(i?)[<.A,A/B(i?) + t^A,MBiE)]\MB > , (69) 
T^^NMBiE) =< ^i-UE)\hU..GpNmtpN,MB{E) + 

ipN.MB 

{E)]\MB>, (70) 

(E) i;[J^{E)\hU^^G.N{E)[t aN,MBiE) + t^]^j,^g{E)]\MB > . (71) 
In the above equations, the ttttN scattering wave function is defined by 

< = < ^^m^^NiE) , (72) 

where the scattering operator is defined by 

ni-J},{E) n7rN\[l + t^^N,..N{E)—— r^— -] • (73) 

Here the three-body scattering amplitude tTrTTN,TVTTN{E) is determined by the non- 
resonant interactions VTr-n-, Vt^n,-kN and VTrTTN,-K-KN , as defined by Ea.l|66)). 
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Figure 9. Graphical representations of T^^jv MB defined by Eqs. II67I I-I I73II . 

We note here that the direct production amphtude T^^"^ (E) of Eq. ((68|) is due 
to Vt,t:n,mb interaction, while the other three terms are through the unstable ttA, 
pN , and aN states illustrated in FiglHl Each term has the contributions from the 
non-resonant amplitude tM' b' ,mb{E) and resonant term tfj,g, Msi^)- 



5. Cross Sections and N-N* Transition Form Factors 



In this section, we give formula for calculating the cross sections of all electrowcak pion 
production reactions. Their relations with the commonly used CGLN and mutipole 
amplitudes are given in appendix B. For later discussions in section 5, we also present 
formula for calculating the electromagnetic N-N* transition form factors which are 
the main focus of recent studies of electromagnetic meson production reactions. 



5.1. Cross Section Formula 



With the relation Eq. (l46|) between the S- and T- matrices and the normalization 
< k\k >= S{k — k ), the amplitude T^N,TrN for the pion photoproduction reaction 
7(g*) + N{-q) 7r(fc) + N{-k) defined by Eq. gS]) can be written in the final ttN 
center of mass frame as (suppressing spin-isospin indices) 

rr 1 ™^ r 7 1 t'7A\ 

(27r)3 ^ENiq)ENik)2E^ik)2q' 

Here is the polarization vector of photon. In the tree-diagram approximation, 
the current matrix element Jem is of the form of [u-p I U-q\, where I is the usual 
invariant amplitudes calculated from the Lagrangian L{x) ~ where 
Jem (•^) the electromagnetic current operator and (x) is the electromagnetic field. 
Similarly the amplitudes for the electroweak pion production reactions e{pe) + N (p) 
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e'{p,,)+n{k)+N{p'), vM+N{p) -> (pe-) + Ti{k) + N{p'), and v{p,)+N{p) 
^{Pv') + ^{k} + N{p') can be written as 

1 TTlN 



(27r)9/2 ^EN{p)EN{p')2E^{k) q 

1 TO AT Gf COS 



(27r)9/2 ^EN{p)EN{p')2EAk) V2 



Jem ' ^emi ('^^) 
:c • ice, (76) 



rp _ ^ Gf , . 

(27r)9/2 ^ENip)ENip')2E^{k) ^ 

where J^^,, J^^ are the matrix elements of charged current and neutral current, 
respectively. The lepton current matrix elements are 

LJf„=S(pe')7^"(Pe), (78) 

LJf, =S(pe')7^(l-75)«(p.), (79) 

LJ:, =S(p,07^(l-75)^(p.). (80) 

The pion production current, J^(a = em, cc, nc) can be written in terms of commonly 
used CGLN amplitudes and multipole amplitudes. These are summarized in appendix 
B. 

The differential cross sections of pion productions reactions due to electromagnetic 
(em) and charged weak current (cc) in the massless leptons (toe = 0) limit can be 
written as 

c^o-L 1 e-* E^. Q2 ^ "^JV a2 p N 

^ yir~E-l ^ern, (81) 



dEe'dne'dni 4Q^ Ee 1-e 2TT^mN MttE' 



(l^)^i?cc. (82) 



dEe'dQe'dni 2 E„l-e2T:^mN '^■nE 

where E is the invariant mass of the final nN state, e is defined by the lepton scattering 
angle Oi^p as e = 1/[1 + 2 ^"^J tan'^ ^^nd /c^ is pion momentum in the TriV center 
of mass system. The functions Ra depends on the pion angle with respect to the 
direction of momentum transfer q and also the angle between the the tt — N plane 
and the plane of the incoming and outgoing leptons. Explicitly, we have 



i?em = Rln + ^Rem + ^2^(1 + e)i?em,c COS (j)^ + ei?fj,, COS 2(j)^], (83) 

i?cc = + <^Rcc + V2e(l + e)(^cc^c COS (/i^ + R^^^ sin 0^) 

+ e(i?Jj^^ cos 2(f)^ + R^^ sin 2</>^) . (84) 

The structure functions R^ in the above equations are calculated from the current 
J/^ for the N + ja ^ TT + N introduced in Eqs. (l74p - ([77)) in the pion- nucleon center 
of mass system. It is common to choose the momentum transfer of leptons as the 
quantization z-direction q = \q |(0,0, 1) and set the outgoing pion on the x-z plane 
fc^ — Ifc^rKsinfl, 0, cos^). The structure functions can then be written as 

\ TX\2 , \ Ty\2 

Rl - Et °' ^' °' -v/r^im(j-jr)], (85) 
-E^|'^°^ (86) 




RZ = H jS[-Re(J^ JD + \/^Im(J^ JD] , (87) 
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2 




1 -e. 



_ ^^^^_:[Re(J° JD + \/ ^Im( J°Jr)] 

(91) 

where a = em, cc, and we have defined 

Jl = Jl+ojq- J^IQ" . (92) 
The spin sum of the nucleons ^ is 

For investigating the weak pion production reactions induced by [i neutrinos, the 
above formula need to be modified to include the finite mass to^ of the outgoing 
lepton. These formula were given in Ref. [6] and were used in obtaining the results to 
be reviewed in section 6.2. The cross section formula for the neutral current reactions 
can be obtained by replacing Gp cos 9c and Jcc of Eq. ([82| with Gp and Jnc- 

For the structure functions of the electromagnetic current Rem, we use J" = J° 
and Im(Jg^J|'^) = Im(Jg„JJ'^) = 0. For pion electroproduction cross sections, it is 
convenient to write Eqs.([8T|l as 



with 



^ = Tr— (94) 

dEe'dfle'dni dVLl' ^ ' 

where — (E^ — m|,)/(2£') and q-yx ~ {E'^ ~ to^)/(2tojv)- 



5.2. N* Transition Form Factor 

The main objective of analyzing the data of electromagnetic meson production 
reactions is to extract the 'jN N*{JT) transition form factors with J and T 
denoting the spin and isospin of a nucleon resonance. In this section, we define these 
quantities within our formulation. 

Our starting point is the following Lagrangian density within the framework of 
the relativistic quantum field theory 

Lemix) = ej^j^{x)Af^{x) , 
where A^(x) is the electromagnetic field and jemi^) is the current operator. In the rest 
frame of N* , the electromagnetic 'yN{sz,tz) N*{JT) transition form factors are 
usually characterized [44l [45] by the helicity amplitudes A\ for the spatial components 
and 5*1/2 for the time component of currents : 

Mj2,tAQ^) N*{JT)\jem{Q^) ■ ei\N{sz = 1/2, t,) > , (97) 

Mf2,tAQ^) N*{JT)\jern{Q^) ■ ei\N{sz ^ -1/2, t,) > , (98) 

^V2,t.(Q') N*{JTM^iQ^)\N{sz = 1/2, i,) > , (99) 
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where — —q^ = cp — uP' is defined by the photon momentum — {uj,q), and 

The efi^ective photon energy is determined by the resonance mass Mres as — 
(M^gj - m%)/{2Mres)- Thejiehcity ampUtudcs Eqs. (HHl-lIMl) are related to the 
radiative decay width of the N* as 

[Width],,, (TV* (JT)) = ^^^^[\AiJ,,f + \A(J,,f] .(102) 



Since the nucleon resonances couple with the meson-baryon continuum states, the 
N* state vector appearing in Eqs. (|57|) - ([M)) is an eigenstate (Gamow state) of the 
Hamiltonian at the resonance energy Eres = [Mres, — ^r^es/S) which is defined by the 
condition Eres — + Y.{Eres)- It consists of a bare N* state and meson-baryon 

components 

\N*{JT) > = \N*{JT) > 

+ ^ {SmB.M'B' + tMB,M'B'GM'B')^N*^M'B'\N*{JT) > 
MB,M'B' 

= \N*{JT) > +^ >< MB\Tn'-^mb\N*{JT) > . (103) 

MB 

Here we have used the relation Eg. ([55]) for defining the dressed vertex Tn-^mb- 
Thus the form factors defined by Eqs.([97 |) - ((99l) are determined by the following matrix 
elements 

< N*{JTMJN > -e^ =< N*{JTMJN > , (104) 

where the meson cloud effects are 

Smc^Y.< N*iJT)\rN'^MB\MB > Gmb[< MB\j^JN > -e^] . (105) 

MB 

The matrix element [< MB\j^.^-^\N > -e^] defines the non-resonant vmb,^n parts of 
the interaction V22 of Eq. (gS]). Eq. (fTM)) is illustrated in Fig. ^ 

Our normalization is chosen such that the vertex functions F^at^tv* and T^m^n* 
of Eqs. ([5^ - ([55|) in each partial wave are related to the matrix element of the current 
operator by 



< N*{JT)\ 

^3 em 



I2J + 




An 


I2J + 


h 



< N*{JT)\ejem • > - y ^^F^jv^at. ( JT) . 

For comparing with theoretical predictions from hadron models and LQCD, we 
need to evaluate the helicity amplitudes Eqs. (|97p - ([M|) at the resonance pole E^es- 
This is a non-trivial problem and is being investigated in Ref. . 
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Figure 10. 
Eos.l fTOit -l lTOSl l 



Graphical representation of tiie dressed F-yjv— »jv* defined by 



6. Results 

With the formulation presented in the above two sections, very extensive data of 
ttN, ^N, N{e,e') and also N{v^,^Tr) reactions have been analyzed. Most detailed 
resultsgl [51 E] are for the A(1232) state. These will be reviewed in subsection 
6.1 for the electromagnetic 'jN — > ttN and iV(e,e'7r) processes and 6.2 for the weak 
N{Vfj_, , /i7r) reactions. The investigation of higher mass N* states began in 2006 and is 
still in the progressing stage. Thus only limited results will be reviewed in subsection 



6.1. Electromagnetic Excitation of the A(1232) state 

The electromagnetic excitation of the A(1232) state was studied in Refs. [ll[5l[9]. 
The main objective was to extract the A(1232) form factors from the data 

of photoproduction and electroproduction of tt in the invariant mass W < 1.3 GeV 
region where only nN and 7iV channels are open. Thus it was studied by using the 
formula presented in section 4 by keeping only one bare A state and including only 
the ttN and 7iV channels. The resulting model is identical to the model developed in 
Refs.U (called the Sato-Lee (SL) model in the literatures). 

The 7iV A (1232) form factor rA,7Ar is parametrized in the form developed 
by Jones and Scadron 46J. With the normalization < k\k' >— S{k — k') for the plane 
wave states and < (t>B\<i>B' >= Sb,B' ior B = N and bare A states, the covariant 
form of Jones and Scadron can be cast, in the rest frame of the A and for the photon 
momentum q = {Lu,q), as 

< mj^,mt^\TA,fN{q)\>^f>^N,mt,~, > 



= F X (|TOt^|ilTOt„0) 

X [M„^.^ ,A,,A„ iq)GM{Q^) + i?™,^ A„ iq)GE{Q^) + C^,^ ,A,A„ {q)Gc (Q')] , (106) 



where < jm\ji, j2,mi, m2 > is the Clebsch-Gordon coefficient of ji + j2 = j coupling, 
A-y and Xn are the helicities of the initial photon and nucleon, ruj^ is the z-component 
of the A spin, nit^ and mt„ denote the isospin components. In Eq. ()106p we have 
defined 



6.3. 




(107) 
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Figure 11. Phase shifts of ttN elastic scattering up to Tj^ =250 MeV. SoUd 
and dotted stand for model SL and SL2 respectively. Data, L2T 2ji are from the 
energy independent SAID |47| analysis plus 8 points from their energy dependent 
solution for the P13 and P31 partial waves at lower energies. 



and the excitation kinematics are contained in 

Mmj^,\-,XNi<l) = < mj^\iS X q- ex^\XN > , (108) 

Emj^,\^\i^iq) = < mj^lS" • e*A^(T • (f + 5" • (f<T • eA^IAjv > , (109) 

Cmj. ,A^A«(g) = — < nijJS ■ qa ■ qcqIXn > , (HO) 



where e = photon polarization vector is defined by e*±i — ^^{x ± iy), 

and £j.i — for = ±1, eq = and £9 = 1 for the scalar component X^ = 0. The 
transition spin S is defined by < jA'mAlSmljN'm'N >—< jA'mAljN^'mN'm >■ 

The form factors Gm{Q^), Ge{Q^), and Gc{Q^) in Eq. (|106p describe magnetic 
Ml, Electric E2, and Coulomb C2 transitions. Choosing the photon direction q in 
the z-direction, the above form factors are related to the form factors in helicity 
representation defined in Eqs.(f97 |l - ([99|) . which are consistent with the convention of 
Particle Data Group [1] (PDG) 

^3/2(Q') = - ^[Gm{Q^) + Ge{Q^)] , (111) 
^1/2(Q') - - |[Gm(Q') - 3G£;(Q2)] , (112) 

5i/2(g2) = - 4|^Gc(Q'), (113) 



with 



Dynamical Models of the Excitations of Nucleon Resonances 



20 



The dressed form factor T/^ -yPf has the same symmetry property of the bare 
vertex defined above. Thus it can be expanded in the same form of Eq. (|106p . We 
denote the dressed form factors by Gm(Q^), Ge{Q^), GciQ"^)- The corresponding 
hehcity amphtudes A\ can also be calculated by using the same relations Eqs. (|llip - 
(|113p . In Ref. [4], it was shown that rA,7Af can also be calculated from the K-matrix 
form of rA,7Af which is directly related to the imaginary parts of the full multipole 
amplitudes Mi+, Ei+ and Si+ at the resonance energy Wr where the TriV phase shift 
is 90'^, independent of the form of the non-resonant amplitudes. Thus the dressed 
ratios can be calculated from 

REM{W = Wn)=^^^'f§^, (115) 

RsMiW^Wn)^^^^^, (116) 
2mA Gm lmMi+ 

It is common to define G%j for the Ml transition form factor which is related to our 
dressed form factor by 

V ^/l + QyimA+mNr 

where = 115 MeV is used in extracting the data from Alf^ amplitude of pion 
electroproduction amplitude and T^^ — 93 MeV from the constructed model. 

With the above definitions of 'yN A (1232) form factors, we now describe the 
results obtained in Refs.jU [5l [9]. The first step in extracting the ■yN A (1232) 
form factors is to fix the hadronic parameters by fitting the ttN elastic scattering up 
to W = 1.3 GeV. Two fits from Refs. 01121 are shown in Fig[Tl] These two models 
will be called SL and SL2 models in later discussions. Their differences are mainly in 
fitting the weak P13 partial waves. These two fits provide us with an opportunity to 
examine the model dependence of the extracted •yN A (1232) form factors. 

The next step is to adjust the bare A (1232) form factors Gm{Q^), 

Ge{Q^), and Gc{Q^) to fit the world data of 7p Tr°p,n~^n , p{e,e'TT°)p and 
p{e, e''K^)n. In Fig ll2[ we show some typical fits to the structure functions of 
p{e, e'n^)p . The resulting bare (solid triangles) and dressed (solid squares) form 
factors are shown in FigfTSl In the same figure we also show the LQCD results 
(open crosses with errors) which are obtained from applying a chiral extrapolation 
procedure to get results in the physical region from the calculations with very large 
quark masses. We see that LQCD results agree only very qualitatively with either the 
extracted dressed or bare form factors. There are several difficulties in interpreting 
these results, as discussed by Pascalutsa and Vanderhaeghen [57]. First, the chiral 
extrapolation is only valid for low Q^, although it has been used in a rather high 
region. Second, there are higher order corrections on the commonly used chiral 
extrapolation, which have not been under control. Thus it is not clear what to conclude 
from Fig. [13] for the results from LQCD of Ref. [55l[56]. Further investigations are 
clearly needed. 

Here we note that the extracted bare form factor Gm(Q^) (solid triangles) in 
Fig. [13] are close to the following parametrization of Ref. [5] 

Gm(Q') = GM(0)i?sL(g')Gp(Q2) , (118) 
where Gp{Q^) = 1/{1+Q'^/Mlf with = 0.71 {GcV/cf being the well determined 
nucleon form factor, and 

i?SL(Q') = (l + aQ')exp(-6Q2)^ (Hg) 
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p(e.e'p)7i° W=1221 MeV Q'=0.06 GeV' 




-G,6 -0-3 -0-e -0.3 -0,6 -1 -0-3 -0-6 

COS i},* COS 1*,* COS iJ,* COS it: 



p{e,e''p)7T° W=1737MeV Q'=0.177GeV' 




COS COS COS iJ,' VV (GEV) 



p(e,e'p)7T° W=1?40 MeV Q'=0.9 GeV^ 




COS COS COS -s; 



p(e.e'p)7T° W=1240 MeV Q'='-45GeV' 




:0S iS.,* COS !>„• COS iS,* 



Figure 12. Fits to experimental p(e,e'p)n^^ structure functions. Solid lines are 
from the fits with the bare form factors GMiQ'^)t Ge{Q'^) and Gc{Q^) adjusted 
at each . The dashed curves are from the calculations using the parametrization 
Eqs.| |118I I-| |119I I. The structure functions aa are -R"^ defined in Ea.l l83l l. Data 
are from MAMI [48] at = 0.O6 GeV^, BATES gOlElEI] at = 0.127 GeV^, 
CLAS [52] = 1220 MeV) and MAMI [53] = 1221 MeV) at = 0.2 GeV^ 
and CLAS [54] at = o.9, 1.45 GeV^. 
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Figure 13. The extracted ■yN —> A form factors. Dark squares (triangles) are tlie 
dressed (bare) values. Open crosses witii errors are the lattice QCD calculation 
of Ref. [g5l[56] . 







Rem{%) 






Rsm{%) 






UIM 


SL 


SL2 


UIM 


SL 


SL2 


0.16 


-1.94(0.13) 


-2.45(0.2) 


-2.57(0.2) 


-4.64(0.19) 


-4.44(0.35) 


-4.36(0.35) 


0.20 


-1.68(0.18) 


-2.21(0.2) 


-2.31(0.2) 


-4.62(0.18) 


-4.23(0.35) 


-4.14(0.35) 


0.24 


-2.14(0.14) 


-2.70(0.2) 


-2.76(0.2) 


-4.60(0.28) 


-4.32(0.35) 


-4.21(0.35) 


0.28 


-1.69(0.27) 


-1.99(0.2) 


-2.07(0.2) 


-5.50(0.31) 


-5.08(0.35) 


-4.97(0.35) 


0.32 


-1.59(0.17) 


-2.29(0.2) 


-2.35(0.2) 


-5.71(0.33) 


-4.87(0.35) 


-4.75(0.35) 


0.36 


-1.52(0.27) 


-1.80(0.2) 


-1.82(0.2) 


-5.79(0.43) 


-4.76(0.35) 


-4.56(0.35) 



Table 1. Extracted values of E2/M1 ratio Rem and C2/M1 ratio Rsm = 
5i+/Mi+ at Q2 = 0.16 - 0.36 GeV^ from analysis of results from a CLAS 
measurement I52| of the p{e, e'p)n'^ reaction. Methods used are Unitary Isobar 
Model (UIM) and the SL and SL2 which use hadronic parameters determined in 
Ref. ^ and Ref. |9], respectively. Errors are statistical only. 



with Gm(0) = 1.85, a = 0.154 (GeV)-^ and b ^ 0.166 (GcV)-^. By using this 
parametrization, the predicted bare (dotted curve) and dressed (sohd curve) G\j{Q'^) 
( defined by Eq. (|117p ) are compared with the available empirical values in FiglUl It 
is clear that the resulting dressed G\.j{Q^) (solid curve) agree well with the available 
empirical values. The differences between the solid and dotted curves indicate that 
the meson cloud effects, illustrated in Fig llOl are important in the low region and 
gradually diminish as increases. This result is one of the main accomplishments 
of many-year study of A^-A (1232) excitation, and has motivated future studies up to 
Q2 = 11 (GeV)2 with 12 GeV upgrade of CEBAF at JLab. 
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(GeV^) 



Figure 14. Magnetic dipole transition form factor G*j^j for 7* A' A(1232), 
normalized to the proton dipole form factor GoiQ^) = + Q'^/A'^]^ with 

= 0.71 (GeV/c)^. Experimental points are analyses of inclusive data (Q) from 
pre-1990 experiments at DESY and SLAG |58l 159 60 61] and recent exclusive 
p{e,e'p)-K° data (blacksquare) from BATES [gj [501 15ll, MAMI [48| |53] and 
JLAB [12 m |62] |63] Ell m]. Solid curve is from the dressed calculation of 
SL model using the parametrization of Eq. I I119I I. The dotted curve is obtained 
when the meson cloud effect, defined by Eo. lllOSII is turned off. 



Historically, the A (1232) is described by the constituent quark model. To see 
the extent to which the extracted Gm{Q^) form factors can be understood with this 
model, it is instructive to first consider the naive s-wave non-relativistic quark model 
within which ^Lp for the proton magnetic moment and /iA+p for the IS.'^-p Ml transition 
are defined by 



2m, 



2m, 



= (A+,m, 



1 1 \ " 

2'^ 2m, 



(Ti(z)\p,nis 



11 V — 

2 ' ^ 2m, 



ai{z)\p,m,s 



(120) 
(121) 



From the above relation and the definition Eq. (|106p . one observes that the magnetic 
Ml form factor of — > A at — q can be directly calculated from the proton 
magnetic moment 



Gm(0) = [\/2Gp(0)] 
where q = {m\ 



2{EN{q)+mN) 



3(mA + mAr) 



0.84/ip . 



(122) 



lr)/2mA ' 

exp 



I ENiq) 

260 MeV/c. If we use the empirical value of proton 



magnetic moment [fp'^ = 1 + Kp ~ 2.77, we then find Gm(0) 2.32 which is 

considerably smaller than the extracted dressed value ~ 3.2 seen in Fig.[T31 This was 
observed in Ref.[lj and interpreted as due to the large meson cloud effects which are 
the difference between the solid and dotted curves in FigfMl 
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We thus observe that extracted bare value Gm(0) = 1.85 can perhaps be 
understood in terms of constituent quark degrees of freedom if we tune properly 
the constituent quark model calculations. On the other hand, our extracted bare 
E2 transition form factor Ge{Q) cannot be understood within the non-relativistic 
constituent quark model. With the tensor force within the conventional one-gluon- 
exchange, the estimated E2 transition of — > A is known to be negligibly small 
compared with the value calculated from our value G'_b(0) = —0.025. In Ref.[9J, the 
extracted form factors are also compared with relativistic constituent quark models. 
Only qualitative agreement is obtained. 

We next present our determined dressed Rem and Rsm in the low region 
where very large meson cloud effects have been identified in Fig. [131 Our results, 
SL and SL2, are listed in table [1] and compared with the values determined using 
the unitary isobar model (UIM). The differences between our values and that from 
the UIM reflect some model-dependence in the extraction. Here we note that only 
the data of five of the eleven A^(e, e'7r)iV independent observables were available and 
used in the fits. Thus the differences between different models shown in Table [T] are 
surprisingly small. So far there is no satisfactory theoretical understanding of the 
results of Rem and Rsm shown in Table 1. 



6. 2. Weak excitation of the A state 

The model developed in Refs.[4l [5], the SL model, was extended to investigate 
neutrino-induced pion production reactions. The extension is tedious but 
straightforward, as detailed in Refs.^iTj. Here we just focus on the extraction of the 
weak iV-A (1232) form factor which has vector {V) and axial vector [A) components. 
The vector current matrix element < A I I N > can be obtained from the SL 
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Figure 15. Total cross sections of N{vn, ^~ ■k)N reactions predicted by the SL 
Model|6j. The data are from Ref. 71.. The soHd curves are from full calculations. 
The dotted curves are from turning off pion cloud effects on N-A transitions. The 
dashed curves are the contributions from the non-resonant amplitude. 
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model by appropriate isospin rotations. The most general form of the axial vector 
current matrix element is well known, as given in Refs.[5ni [SZl EH]. To see how it is 
different from the electromagnetic excitation given in Eqs. (|106p - (|109p . we cast[n| it in 
the rest frame of a A on the resonance energy ( pA = (tia, 0)tPn = [ENiq), —Q), 1 = 
(mA - ENiq),q}) as 



<A\A^\N> ^^lEE±J^[(^d, + ^I^^^d,)S 
" 2mAr mf^ 



(d. + d3)^^-.d44;i;^^]r, (123) 
m^ mj^{EN + 'mN) 



<A\A"'^\N> ^ E^ + mN^^^S-qjm^ + E^) 
V 2mAr mff 

_^^S^q{rn^-E^^^,^ (124) 
m^ 

where T* is the z— th component of the isospin transition operator (defined by the 
reduced matrix element < 3/2 || f \\ 1/2 - < 1/2 || f+ || 3/2 >= 2 in Edmonds 
convention [69] ) . and the transition spin S is defined by the same reduced matrix 
elements of T. The above expression suggests that di , d2 terms describe the Gamow- 
Teller transition and d/i describes the quadrupole transition. For simplicity, we follow 
Ref. |68j to fix the form factors di{q^) a,t — using the non-relativistic constituent 
quark model. The axial vector current operator for a constituent quark is derived 
from taking the non-relativistic limit of the standard form gAqqi^lb'^q- By some 
derivations [B], we find that 

d,{Ql) = 9\{Ql){l + , "f-"^^ ) , (125) 
d2{Ql) = - gl(Qo) ^. , , (126) 

2 

dsiQl)^ -9*AiQl)^^, (127) 

where g^iQa) = ^Tf t^^ with gA = 1-26 and Qq = (mA — m^r)^. This agrees with the 
results of Ref. ^68] if we neglect the difference between mjv and mA. 
To account for the g^-dependence, we assume that 

d,{Q^) = d,{0)RsL{Q^)GA{Q^) , (128) 

where Rsl{Q^) is defined in Ea. (|119p and has been determined in the study of 
jN ^ A (1232) form factor, and Ga(Q^) = 1/(1 + Q^/Ml)^ with Ma = 1.02 GeV is 
the nucleon axial form fact or [70]. 

With the axial form factors defined above, our calculations of p(^'^t, /i7r)A^ do 
not involve any adjustment of the parameters, since all of the the parameters of the 
non-resonant amplitudes and the vector part of the A^-A transition form factor have 
been completely fixed in the study of electromagnetic pion production. The predicted 
total cross sections are compared with with the data[7T] in FigdS] We see that the 
predictions (solid curves) agree reasonably well with the data for three pion channels. 
For the data on neutron target, our predictions (solid curves in the middle and lower 
figures) are in general lower than the data. This is perhaps related to the procedures 
used in Ref. [71] to extract these data from the experiments on deuteron target. 
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Similar to the electromagnetic A^-A transition, we have also found significant 
meson cloud effects on the axial A^-A transition form factor. This is also shown in 
FigdU We see that our full calculations(solid curves) are reduced significantly to 
dotted curves if we turn off the dynamical pion cloud effects. If we further turn 
off the contributions from bare N-A transitions, we obtain the dashed curves which 
correspond to the contributions from the non-resonant amplitudes. Clearly, the non- 
resonant amplitudes are weaker, but are also essential in getting the good agreement 
with the data since they can interfere with the resonant amplitudes. 
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Figure 16. Differential cross sections da/dQ^ of p{u ^ , ^~ ■K'^)p reaction averaged 
over neutrino energies 0.5 GeV < Ev < 6 GeV. The curves are the predictions 
of the SL Model[B]. The dotted curve{dot-dashed curve) is the contribution 
from axial vector current A (vector current V). The solid curve is from the full 
calculations with V-A current. The data are from Ref. |71). 

In FigHniwe compare the Q^-dependence of the differential cross sections da/dQ^ 
with the data from ANL[7T]. We see that our predictions(solid curve) agree reasonable 
well with the data both in magnitude and Q^— dependence. In Fig[Tn] we also 
compare the contributions from vector current (dot-dashed curve) and axial vector 
current(dotted curve). They have rather different Q^-dependence in the low region 
and interfere constructively with each other to yield the solid curve of the full results. 
Since vector current contributions are very much constrained by the (e, c'tt) data, the 
results of Fig [11] suggest that the constructed axial vector currents are consistent with 
the data. 

The extraction of the axial iV-A form factor is much more difficult because the 
lack of sufficient data. The dressed (solid curve) and bare (dotted curve) axial A^-A 
form factors are shown in the right-hand side of Fig ll7l Clearly, their Q^-dependence 
is weaker than the — > A form factors which are discussed in the previous subsection 
and also displayed in left hand side of FiglTT] However, the meson cloud effects, the 
difference between the solid and dotted curves, are comparable in both form factors. 

The axial A^-A form factor was determined in previous analysis. In Fig[T8l we 
see that our results (solid) are significantly different from the previous results ( dot- 



Tot 
V 
A 




Dynamical Models of the Excitations of Nucleon Resonances 



27 



dashed curve) at high Q^. Obviously, more experimental data are needed to resolve 
the differences. With the new world-wide effort in developing next-generation neutrino 
experiments, progress in this direction is expected in the near future. 
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Figure 17. The N-A form factors: left panel: Magnetic Ml form factors given 
in Ref. 5 , right panel: axial vector form factor determined in Ref. [6]. The solid 
curves are from full calculations. The dotted curves are obtained from turning off 
the pion cloud effects. Go = 1/(1 + Q^/M^)^ with My = 0.84 GeV is the usual 



proton dipole form factor and Ga = 
the axial nucleon form factor of Ref. 



1/(1 + Q /M\Y with Ma = 1.02 GeV is 
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Figure 18. Compare the dressed axial N-A form factor predicted by the Model 
of Ref. j6j (solid curve) with the empirical form factor (dot-dash curve) determined 
in Ref. [72). 



6.3. Excitations of higher mass N* states 

To investigate higher mass N* states up to invariant mass W = 2 GeV, we apply the 
full model developed in sections 3 and 4. The meson-baryon (MB) channels considered 
are 7A'^, ttA'^, 'qN and the ttttA^ channel which has resonant ttA, triV, pN components. 
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Figure 19. The ttN partial wave amplitudes of isospin T = 1/2 calculated 
from the JLMS model |TOl are compared with the energy independent solutions of 
Ref. [47) . Left (right) panel is for real (imaginary) parts of the amplitudes 



The resonant amplitude i^/^/ mb Eq.igS]) are generated by including one or two 
bare N* states in each partial waves. Clearly it is a highly nontrivial task to extract 
the resonance parameters from solving this multi-channels multi-resonance problem. 
It requires simultaneous fits to all available data of ttN , jN, and A^(e, e') data with all 
possible two-particle and three-particle ttttN states. This ambitious work started in 
2006 at the Excited Baryon Analysis Center (EBAC) of JLab, and is still progressing 
rapidly. Thus the results reviewed in this subsection are only the first-step results 
which will be refined when all of the world's meson production data of TriV, 7iV, and 
N{e, e') reactions are included in the analysis. 

6.3.1. ttN scattering Similar to the study of the A (1232) state, the first step to 
investigate higher mass N* states is to determine the hadronic parameters by fitting 
the data of ttN elastic scattering. Such a fit was obtained in Ref.[Tn] by assuming one 
or two bare N* states in each of S, P, D, and F partial waves . The ttN scattering 
amplitudes of isospin T = 1/2 predicted by the resulting model, the JLMS model, are 
compared with the empirical values of SAID [47] in Fig[T9l Similar good agreement is 
also found for the T = 3/2 partial waves, as also given in Ref. JIOI. The corresponding 
good agreement with the data of differential cross sections and polarization observable 
P are illustrated in FigI5D]for some of the data. The predicted total cross sections are 
also in good agreement with the data as shown in Fig|2T] 

The resulting parameters of 21 bare N* states, presented in Ref. [TU] , is the starting 
point for performing a dynamical coupled-channel analysis of the world's meson 
production data of ttN, 'yN, and N{e, e') reactions. In the next three subsections, we 
review the results obtained so far. Here we also mention that it is necessary to develop 
an analytic continuation method to identify the nucleon resonances with the poles of 
the scattering amplitudes on complex energy plane. This has been developed [T4]. but 
will not be discussed here because of its technical complexities. 

6.3.2. ttN ttttN reactions The main difficulty in fitting the ttN elastic 
scattering data, described above, is that the model contains many parameters 
mainly due to the lack of sound theoretical guidance in parametrizing the bare 
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Figure 20. DifTerontial cross sections dcr/dQ (left) and asymmetry P (riglit) of 
TT^p n^p,n^'n reactions. The solid curves are from JLMS model [To], 

N* ttN, r/N, ttA, pN, aN form factors. Thus it is necessary to examine these N* 
parameters; in particular the parameters associated with the unstable ttA, pN, and 
aN channels. This has been done in Ref.fTT in the study of ttN ttttN reactions 
which arc known to be dominated by these unstable particle channels. 




' W (GeV) ' ' W (GeV) 



Figure 21. Total cross sections of Tr+p (left) and tt~p (right) reactions. Solid 
curves are from the JLMS model [lO]. Only few data are shown for a clear 
comparison between curves. The data are from Refs. [Tl 1731 . 

Before we present the predicted TriV — s- ttttN cross sections, we note here that 
the main feature of our approach is a dynamical coupled-channels treatment of the 
unstable TrA,pN,aN channels. This effect can be explicitly seen by writing the 
coupled-channels equations, Eq. (|57|) . as 

tMB..NiE) = E [1 - ^^]Mi3.M's'^M'S',.7V , (129) 

m' b' 

where MB = ttA, pN, aN, and the intermediate meson-baryon states can he M B = 
TrN,r]N,TrA,aN, pN. The predicted TriV wttN total cross sections depend on the 
coupled-channel effects due to these intermediate M'B' states. 

The results for irN ttttN total cross sections are shown in Figl22l We see that 
our full calculations (solid curves) can reproduce the data to a very large extent for all 
possible ttttN final states up to = 2 GeV. These results are far more successful than 
all of the previous investigations, as discussed in Ref.fTSj. When only the term with 
M'B' = MB in the Eq. (fT29l) and in Tn'^mb of Eqs.(l52|-(l53l) is kept, the calculated 
total cross sections (solid curves) are changed to the dotted curves in Figl22l If we 



Dynamical Models of the Excitations of Nucleon Resonances 



30 




Figure 22. The coupled-channels effects on vrA'^ — > ■kttN reactions. The 
solid curves are from full calculations, the dotted curves are from keeping only 
M'B' = MB'm the Ea.| [T29l l and inV m-'-.mb of Eqs.l|52ll-((53j, the dashed curves 
are from setting t},iB.M' B' — ^MB.M' B' ■ The data are from 1741 . 



further neglect the coupled-channels effects by setting t-^NMB = V7rN,MB, we then get 
the dashed curves which are very different from the full calculations (solid curves), in 
particular in the high W region. Clearly coupled-channel effects are very large. 

The results shown in Figl^H indicate that the N* ttA, pN, aN determined from 
fitting ttN elastic scattering data are reasonable, but clearly need to be improved. To 
make the progress in this direction, it is necessary to have more complete data of 
ttN — > ttttN reactions from new hadron facilities such as J-PARC. Hopefully, this can 
be realized in the near future. At the present time, we have to rely on recent data of 
7iV ttttN to refine the N* — > ttA, pN,(TN parameters. Effort in this direction is 
being made at EBAC. 

6.3.3. Electromagnetic pion production reactions The fits to t:N reaction data, 
presented in the previous two subsections, have fixed all of the hadronic parameters 
of the effective Hamiltonian Eqs. (l551) - (|15)l . Most of the electromagetic parameters 
associated with the nonresonant — > nN are also known from previous investigation 
of A (1232) state. Thus the bare helicity amplitudes, A3/2, ^1/2, and 5*1/2, 
defined in Eqs. (|97p - (|99|1 . are the main unknown parameters in our investigations of 
electromagnetic pion production reactions. The first step in determining these helicity 
amplitudes had been completed in Ref.[TT] by performing x^— fits to the available 
photoproduction data of 'jN ttN reactions up to = 1.65 GeV. The quality of 
the resulting fit can be seen in Figs. [23] for the 7p n^p. Similar good agreement 
was also obtained for the 77? — > 7r+n, as also presented in Ref. 11 . 

Clearly, the fit to the data needs to be improved, but is sufficient for revealing 
the coupled-channels effects in a dynamical approach. In electromagnetic pion 
productions, the coupled-channel effects are in the loop integrations over the 
intermediate meson-baryon states MB in the following expressions for the non- 
resonant amplitudes and the dressed 7iV N* vertex 

t-KN,-iN = VirN.-yN + ^ tT^N bG mbVM B ,-/N , (130) 
MB 
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Figure 23. Differential cross section da/dQ (upper) and plioton asymmetry 
E-y (lower) for 7p — > n^p calculated from JLMSS model [TT] are compared to 
experimental data obtained from Ref. I47| . 



^N',jN = ^N',iN + ^N*,MbGmbVmB,iN ■ (131) 
MB 

We show the coupled-channels effects on the total cross sections of — > w'^p, 7r+n 
in Fig. [24l We see that the calculated total cross sections (solid curves) are in 
good agreement with the data. The dashed curves are obtained when the channels 
MB = TjN, ttA, pN, and aN are turned off in the loop integrations of Eqs. p30l) - p3Tl) . 
Clearly, the coupled-channels effects 77V rjN, ttA, pN, aN ttN can change the 
cross sections by about 10 - 20 % in the A (1232) region and as much as 50 % in the 
W >1400 MeV second resonance region. 

The meson cloud effects, as illustrated in FigfTUl on several low-lying nucleon 
resonances are also investigated in Ref.|ll|. In general, the resonance parameters 
must be rigorously defined by the poles on the unphysical sheet of complex energy 
plane. This is still being pursued[14j. Here we only illustrate the meson cloud effect on 
the 'yN ttN multipoles for the D13 partial wave. The results are shown in Fig|25l 
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We see that the predicted muhipole ampUtudes agree well with the empirical values 
of SAID [¥7] . and show typical resonant shape at ~ 1.5 GeV. Our model thus also 
has identified a resonance at position close to the Af*(1520, -D13) listed by PDG. If we 
turn off the meson cloud effects on the — > N* in this partial wave, we then get 
the dashed curve. Clearly, meson cloud effects are very large. 

The results reviewed in this subsection are from the very first step of performing 
a dynamical coupled-channel analysis of tt photoproduction and electroproduction 
reactions \xp to W = 2GeV. In parallel, the investigation of ttN ttttN described 
in subsection 5.2 has also been extended to investigate 7*A^ ttttN reactions. Only 
when the world's data of ttTV, 7*7V ttN^ttttN are all included in the analysis, we 
can establish the N* spectrum and their decay properties with confidence. Progress 
in this is being made at EBAC. 




Figure 24. Total cross sections from JLMSS model [TT]. The dashed curves are 
obtained from turning off all MB channels except the ttN channel in the the loop 
integrations in the non-resonant amplitude and the dressed 'yN — > A'^* vertex. 
The dotted curve is obtained by neglecting the off shell effects in the ttN only 
calculation. Experimental data are from Ref. 1471 . 
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Figure 25. The predicted ')N ttN multipole amplitudes in D13 are compared 
with the empirical values of SAID 47 . 
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7. Summary and future developments 

In this article, we have reviewed the dynamical model developed in Refs. [U[5l[6l[71[8l[9l 
[TOl[TT] for investigating the excitations of N* states in irN, -fN and N{e, e'lr) reactions. 
The model Hamiltonian was constructed by using a unitary transformation method, 
and had been used to construct a multi-channels and multi-resonances reaction model. 
The channels considered are 7iV, ttN, rjN, and ttttN which has resonant ttA, pN, and 
aN channels. The resonant amplitudes are generated from 21 bare N* states which 
are renormalized by meson-baryon scattering as required by the unitary condition. 
The model is reduced to the well-studied Sato-Lee (SL) model when only one bare A 
state and irN and jN channels are kept in the formulation. 

The detailed investigations [H [5l [6l [7] of the A (1232) have determined the 
electromagnetic 7iV A (1232) and the axial AN A (1232) form factors. The 
meson cloud effects on these form factors are found to be very large in the low 
region and decreases with Q^. These form factors can be considered along with the 
nucleon form factors as benchmark data for testing the predictions from hadron models 
with effective degrees of freedom and LQCD. 

The investigation of higher mass N* states is based on the full model presented 
in sections 3 and 4. The N* parameters can be reliably determined only when all 
of the available data of ttN, '^N and N{e, e') reactions with all possible two-particle 
and ttttN final states are fitted simultaneously. This ambitious work, started in 2006 
at EBAC, has been progressing well to obtain good fits to the data of ttN elastic 
scattering, ttN ttttN, and 77V ttN reactions. Important coupled-channel effects 
have been revealed. Large meson cloud effects on 7A'' N* have also been identified. 
But more works are needed to establish the extracted N* parameters. 

The current effort at EBAC is to obtain fits to the world data of ttN,j*N 
ttN, rjN, ttttN. Staring with the resulting N* parameters, we then focus on the 
W > 1.7 GeV region by also fitting the world data of nN, -/*N KA, KT,, luN. The 
numerical strategies for handling these additional channels have been developed and 
tested. This effort is needed to face the challenge from the complete and over complete 
measurements of all independent observables of the electromagnetic production of 
KA, KTi reactions. These measurements are expected to be carried out in the next 
few years at JLab. Similar complete experiments are also being developed at Mainz 
and Bonn. 

To end of this article, we point out that the ttN data are very limited except 
the ttN elastic scattering. This could be the main source of the uncertainties of the 
extracted resonance parameters. It will be highly desirable, if more ttN reaction data 
can be obtained at new hadron facility J-PARC in Japan. 
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Appendix A. Interaction Lagrangians 

The expressions of the full Lagrangians for developing the multi-channel multi- 
resonance reaction model are given in Appendix A of Ref. [8]. In this appendix 
we only give the interaction Lagrangians for developing the SL model of electroweak 
pion production reactions. The Lagrangian with tt, p, u, N , and A fields are 

LttNN = - '^''^'^ i'Nlfj.lbfiJN ■ d'^tpT, , (A.l) 
LpNN = gpNN^Nb,j. - '^^'^t^^^'^^P^ ' ■^^^ ' (^-^^ 

ipTTTT = gp7TA4>7r X ^^^Tt] ' , (A. 3) 

Lu;NN = gwNNi^Nbp - TT^^M^'^'^l'^'^V'JV , (A. 4) 

2m AT 



firNA 



LttNA = tp'^TlpN ■ dpCjjjr ■ (A. 5) 

The effective Lagrangians for the Icpton induced electroweak meson production 
reaction are given as 

Leff = -^{-eTe.)Jem,^L 

Gpcos9c,_ ,, N .+ - \ ■ 1 

[^el^{l - ToWcCt. + ^7*^(1 - l5)VeJCcA 

- ^[^>e7^(l - l^>e + e(25^7^ - 25^7'^75)e]jArap. (A.6) 

where a = 1/137, Gp = 1.1664x10-5 GeV-^ and = -1/2 + 2 sin^ 6*^, 9 a = -1/2- 
The Weinberg angle 9\y is known empirically to be sin^ 9\y = 0.231 and cos^c = 0.974 
is the the Cabibbo-Kobayashi-Maskawa (CKM) coefficient. 
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The electromagnetic current (j^,„), weak charge current (j^f,) ^'^'^ weak neutral 
current (j^j,) are written with the iso- vector vector current V^^, axial vector current 
A'^ and iso-scalar vector current V^^^ as 

J'en. = V,^ + Vt, , (A.7) 

J'cc ^{yr + iV^^)-{A't + iA>^), (A.8) 

j,'^, = (1 - 2 sin2 ew)f,,r. - Vt, - . (A.9) 

Here we have neglected the strangeness content of the nucleon. The iso- vector vector 
current and iso-scalar vector current are V^" 

- F2V T -> 

Vfj, = 1pN[FlVlti ~ ^"^"^'2^'^ + 07r X 

+ ^-^[{^n1^ToTxPn) X 0.] + ^e„^^50.(a^^^)a" , (A.IO) 



= ^n[Fis1,. - T^^^^^d^T^^N + ^ea^^s$. ■ . (A.ll) 



F2S 1 , 5pir7 

2m N 2 m. 

The axial vector current needed to construct our model is given as 

1^ = gANl^ls^N - fp^AP^ XTT- F^d^'iT. (A.12) 

Here i^^ = 93 MeV is the pion decay constant, and gA = 1-26 is the nucleon axial 
coupling constant. The iso- vector vector iVA transition current are parametrized in 
the following form 

V,= - i4>''^Tl^f^N + {h.c.) . (A.13) 
The matrix element of iVA current between an N with momentum p and a A with 
momentum can be written explicitly as 



r 



y _ mA -I- ruN 1 



2mN {niA + ttin)^ — q'^ 
12 

{niA ~ niNr - T 

+Gci757 ^-T^ -q^{q^P,-q- Pq,)], (A.14) 

(mA - - T 

The expression for the A'' A transition axial vector current is given in Eqs. (|123|) - 

(UMi). 

Appendix B. Multipole amplitudes of the pseudoscalar meson production 

Here we summarize the formula related the matrix elements (a = em, cc, nc) in 
Eqs. ([55)) - (|92p to the CGLN amplitudes F^ and and multipole amplitudes. Recovering 
the spin indices of J/^, we have 

xlF^Xs = - £1 < -^NisV^Wis) > e„ (B.l) 

where s, s' are spin quantum number of nucleon. F^ is further written as sum of the 
contributions of vector and axial vector currents. 

Fern = fJ^ , (B.2) 

Fee = F^e - F^^ , (B.3) 

Fnc = F^^ — F^^ . (B.4) 
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For electromagnetic reaction, the amplitude(Fem) amplitudes is related to the 
CGLN amplitudeHij as 

The amplitudes are related to those of Ref. [66l as 

= ^F^iAdler), (B.6) 

where E is center of mass energy of pion-nuclcon system. 

The spin structure of the vector E^ and axial vector E^ amplitudes in Eas. (jB.2p - 
(jB.4l) for each of em, cc, nc currents can be parametrized as 

E^ — —id ■ ej_E^ — a ■ ka ■ q x ej_E2 ~ id ■ qk ■ eA_E^ — id ■ kk ■ e±EY 

—id ■ qq ■ eE^ — id ■ kq ■ eE^ + id ■ keoE^ + id ■ qt^E^ , (B.7) 
where e± — q x [e x q) and 

F-^ = -id ■ kd ■ ej_E^ - d -qx e±E.f - id ■ kd ■ qk ■ e^E^ - ik ■ e^E^ 

-id ■ kd ■ qq - eErf - iq ■ eE^ + ie^E^ + id ■ kd ■ qeoE^ . (B.8) 

Here q and k are momentum transfer to nucleon and pion momentum in the center of 
mass system. We defined E^ simply as d ■ kE^ . 

Finally the amplitudes E^ , E,^ are expressed in terms of multipole amplitudes 

TpV,A i,jV,A riV,A tA 



FY 


= + + ipWxMY^ + + ^)p'i-xMY-\ 


, (B.9) 


FY 


= Y,[ii + i)P!MY+ + iP!MY_], 


(B.IO) 


FY 




(B.ll) 


FX 


= J2[-pI'eY+ pI'eY + pI'mY+ - pI'mY^] , 


(B.12) 


FY 


= J2[il + l)PU,Ll-lPY,Ll], 
I 


(B.13) 


FY 


= 5][-a+i)p/Lr++/p/ir_], 


(B.14) 


FY 


^J2[-(i + i)PlsY+ + iP!sY^], 


(B.15) 


FY 


^Y.^{i + i)PU,sY+-iPL,sY,], 


(B.16) 



and 



E,^ = 5][P/< + P!Ef_ + {l + 2]P;M^ + {l- 1)P/M,-1] 
F2^ = + l)P/+iM,t + IPUM,±] , 

= J2[Pl'Ef^ + Pl'Et + Pl'Mfi Pl'Mf'_] , 



B.17) 
B.18) 
B.19) 
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^4 


- Vf-p' 


' pA _ pii 






Fi 




+ ^)P!Lf+ 4 




(B.21) 


Fi 




l)PUiLt - 


-IPl-iLt], 


(B.22) 


F,^ 






-lP!-iSf_], 


(B.23) 


Fi 


-E[-(^ 




ipisn 


(B.24) 



Pl{x) is Legendre function and x — k ■ q. In addition to the normalization of the 
amphtude it is noticed that L^, Sj^^. differ from those of Adlcr. 

The multipole amphtudes are easily calculated from the helicity-LSJ mixed 
representation (Eqs. (C.l) and (C.2) of Ref. :8,:). We express 

<]±\F^\X,Xn> = -^<in/2)j\J^-ex\XN>, (B.25) 

where j± ~ j ± 1/2. The partial wave expansion of the pion production current is 
given as 



< Gl/2)j|Ja-eA|A^ >=2^^y d(cos0)j|^(/,O,l/2,-AW|j,-AW) 

X < 7r(fc), 7V(-fc, s'j^ = -A'^)| J„ • ex\Ni-q, sn = -Xn) > d^^l. (0). (B.26) 



Here we have chosen q = \q\{Q,0,l), k = |A;|(sin 6*, 0, cos 0) and e^j^ = 
(0,Tl/V2, -i/V2,0), e[; = (0,0,0,1) and e'^^ = (1,0,0,0). After some derivation, 
we obtain the following relations: 



Ft = ^-^^[< ..1^-11,1/2 > -Vt^ < ..|F-|l,-l/2 >] , (B.27) 



El = < J-\F'\l, 1/2 >-^J^< J-\F'\l, -1/2 >] , (B.28) 



Ml = ^-(^[< UIF^'I^, 1/2 > +^^^ < j+|P^|l, -1/2 >] , (B.29) 



Ml = J-I^^'ll, 1/2 > -^[^ < j-|F^|l, -1/2 >] , (B.30) 



i^r. = -^-^<. .1^-10, -1/2 >, (B.31) 



ir- = ^<J-l^^''|0,-l/2>, (B.32) 



^r. = ^-^<..|F-|0.,-l/2>, (B.33) 



^r--^<J-l^^^|0*,-l/2>, (B.34) 
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and 



= 4.m{l + l) ^^ ^+1^'''^' > +V ^ < -1/2 >] , (B.35) 



^[-<i_|F^|l,l/2>+Y/[^<i_|F^|l,-l/2>], (B.36) 



^.'t = ^-^<i.|i^^|0.-l/2>, (B.41) 
<i-|F^|0*,-l/2> . (B.42) 



